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Abstract 

A generalized Green's function theory is developed to simulate the inelastic electron tunneling 
spectroscopy (IETS) of molecular junctions. It has been applied to a realistic molecular junction 
with an octanedithiolate embedded between two gold contacts in combination with the hybrid den- 
sity functional theory calculations. The calculated spectra are in excellent agreement with recent 
experimental results. Strong temperature dependence of the experimental IETS spectra is also 
reproduced. It is shown that the IETS is extremely sensitive to the intra-molecular conformation 
and to the molecule-metal contact geometry. 
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Inelastic electron tunneling spectroscopy (IETS) has been recently applied to the molec- 
ular junctions^, 0, [J. The measured spectra show well-resolved vibronic structures corre- 
sponding to certain vibrational normal modes of the molecule. It helps to understand the 
vibronic coupling between the charge carriers and nuclear motion of molecule. The IETS 
has also been found to be sensitive to the experimental setups, for instance, distinct dif- 
ference has been revealed between the spectra of alkanethiol molecular junctions measured 
by Kushmerick et al. and Wang et al. It is therefore expected that the IETS can 
be used as a tool to identify the geometrical structures of the molecule-metal interface, an 
important factor in the design and characterization of molecular devices. The basic theory 
of IETS has been known for many yearsO . Its extension to the molecular wiresj^, Q, atomic 
metal wiresQ] and molecular junctions^, ^, E|, however, has just appeared in the literature 
very recently, partly motived by the success of the experiments. The theoretical modellings 
are generally based on the Green's function theory at different levels, and have contributed 
to many fundamental understandings of the inelastic electron tunneling processes inside a 
molecular /atomic junction. However, none of the theoretical studies has been able to make 
realistic comparisons with the experimental spectra. 

In this letter we present a new computational scheme based on our early developed 
quantum chemical approach for electron transport in molecular j unction (10|. The vibronic 
coupling is introduced by expanding the electronic wavefunction along different vibrational 
normal modes. We have applied this new method to the gold-octanedithiolate-gold junc- 
tion. The electronic structures and vibrational modes are calculated using the hybrid density 
functional theory (DFT) . The simulated spectra are in excellent agreement with the experi- 
mental ones. Our computational scheme can be easily applied to other systems. Theoretical 
simulations are extremely useful for assigning the experimental spectra and to reveal many 
detailed informations that are not accessible in the experiments. 

Our approach is also based on the Green's function formalism. We divide the molecular 
device into three parts, source, drain and extended molecule, as shown in Figure The 
source and drain are described by an effective mass approximation (EMA), while the ex- 
tended molecule is treated with the hybrid density functional theory. The extended molecule 
is in equilibrium with the source and drain through the line up of their effective Fermi level. 
This approach is proved to lead to very good descriptions for the experimental results[10]. 
The current density for the three-dimensional electrodes along the molcular axis z is given 
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where n s and n D are the density of states (DOS) of the source and drain, respectively. Ef is 
the Fermi energy, and T is the device working temperature. Vd is the external voltage. The 
transition matrix element from the source to the drain, T(Vd, Q)T(Vd, Q), is dependent on 
the vibronal motion, Q, and can be written as[10| 

nvo, Q) = E E VMQ)VMQ) E ( J(Q) KiQ>) P) 

where J and K run over all atomic sites of the molecule, which are denoted as 1,2,..., A/", 
site 1 and iV are two end sites of molecule that connect with two electron reservoirs. Vjs 
{Vdk) represents the coupling between atomic site J (K) and reservoirs S (D). Orbital | rj) 
is the eigenstate of the Hamiltonian (Hf) of a finite system that consists of the molecule 
sandwiched between two clusters of metal atoms: Hf \ rf) — e v \ rj). The product of two 
overlap matrix elements (J \ r])(r] | K) represents the derealization of orbital | 77). Here 
parameter z v is a complex variable, z v = E v + iT'L , where E n is the energy related to the 
external bias and Fermi level, at which the scattering process is observed. is the escape 
rate which is determined by the Fermi Golden rule 

rf (Q) = nn s Vj s (Q) \ (J(Q) \ V (Q)) | 2 +7m D V% K (Q) \ (rj(Q) \ K{Q)) | 2 (3) 

All the key parameters are obtained from the calculations of the finite metal-molecule clus- 
ter. The use of the finite system allows us easily including the vibronic coupling since the 
description of vibronic coupling in a molecular system has long been established. The nu- 
clear motion dependent wavefunction can be expanded along the vibrational normal mode 
using a Taylor expansion. The IETS experiment is often done at electronic off-resonant 
region. The adiabatic harmonic approximation can thus be applied, therefore, only the first 
derivative like needs to be considered, where Q a is the vibrational normal mode a of 

the extended molecule. The working formula we have used for IETS calculations follows the 
same principle as the one given by Troisi, Ratner, and Nitzanp. 

We have calculated an octanedithiolate, SCsHiqS, embedded between two gold electrodes 
through S-Au bonds. The extended molecule consists of two gold trimers bonded with an 
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octanedithiolate molecule. The gold trimers were tested in two basic conformations - a 
triangle and a linear chain, see Figure ^ For the triangle configuration, the sulfur atom 
is placed above the the middle of the triangle, resembling the hollow site of a Au(lll) 
surface. Geometry optimization and electronic structure calculations have been carried out 



for the extended molecules of different configurations at hybrid DFT B3LYP level [ll| using 
Gaussian03 program package with the LanL2DZ basis set In the case of the triangle 
configuration, two different schemes have been used for the geometry optimization. The first 
geometry is obtained by optimizing octanedithiol in gas phase and then replace the terminal 
hydrogens with the gold contacts without further optimization (Trl). The second geometry 
is obtained by optimizing the first geometry with the same fixed gold and sulfur distance 
(Tr2), which is also the approach for the chain configuration (Chi). The S-Au distance in 
all calculations is fixed to 2.853 A. The QCME program[13j has been employed for all the 
IETS calculations. 

The calculated IETS of the octanedithiolate junction with the triangle gold trimers are 
shown in Figure El together with the experimental spectrum of Wang et al. at temperature 
4.2 K Q. One can clearly see that the calculated result for Trl configuration is in excellent 
agreement with the experiment. The calculations do not only reproduce all the major 
spectral features observed in the experiment, but also provide very detailed features that 
are smeared out by the background due to the encasing SiaN 4 marked with stars Q|. Both 
theory and experiment shows that the intensity of the vibronic feature follows the order: 
u(C-C) (131mV) > 7(CH 2 ) (155m V) > 5(CH 2 ) (185mV). It should also be mentioned that 
theory and experiment agrees well for the current- volt age characteristics (I-V) which are 
determined mainly by the elastic scattering, see the inserts of the Figure El (A) and (B). 

The IETS calculated with the Tr2 configuration, Figure 121 (C), shows also very rich 
structures. However, its spectral profiles are quite different from those obtained from the 
Trl configuration, as well as the experiment of Wang et a!2\. For instance, the spectral 
peak at 155 mV (7(CH2) mode) has the largest intensity, instead of the mode of u(C-C) 
at 131 mV. The Trl and the Tr2 configurations have the same Au-S bondings, but slightly 
different intra-molecular conformations. The sensitivity of the IETS with respect to the 
molecular geometry is really high. It can thus be concluded that the molecular geometry 
and contact configuration of the device in the experiment of Wang et alQ] is very close to 
our Trl configuration. 



4 



The calculated temperature dependence of the IETS of the Trl configuration is shown 
in FigOJl together with the experimental results of Wang et al|2|. The agreement between 
the theory and the experiment is more than satisfactory. The evolution of the spectral 
bands upon the increase of the temperature is the same for both the experiment and the 
calculation. As an example, in both cases, the peak for mode <5(CH 2 ) at 185 mV disappears 
at 35K, and the peak for mode 7(CH2) at 155 mV becomes invisible at 50K. Our simulations 
indicate that the observed temperature dependence is mainly due to the changes of the Fermi 
distribution. 

In Table HJ the assignments of the vibronic bands observed in theoretical and experimen- 
tal spectra are shown. It can be seen that the vibrational frequencies given by the B3LYP 
calculations are in very good agreement with the experiments. Our computational scheme 
also allows to calculate the spectral linewidth directly, which is determined by the orbital 
characters and molecule-metal bonding, see Eq. |31 The calculated full width at half maxi- 
mum (FWHM) for the spectral profile of mode u(C-C) at 132 mV is found to be around 4.3 
meV, in good agreement with the experimental result of 3.73±0.98 meV^j. However, it is 
also noticed that such a band is attributed from several modes of the same character with 
vibration frequencies 1012, 1038, 1061, 1069, 1084, 1090, 1097, and 1100 cm" 1 , respectively, 
covering a range of 11.0 meV. The overlapping between different vibration modes makes it 
impossible to determine the actual intrinsic lindwidth of the spectral profile from a single 
vibration mode. 

We have calculated the IETS of the gold chain configuration (Chi), shown in FigurelljB), 
to examine the dependence of the IETS on the molecule-metal bonding structure. Indeed, 
the IETS of Chi shows a distinct difference in spectral intensity distribution from those of the 
Trl and the Tr2 configurations. The spectral peak of mode <5(CH 2 ) at 382 mV has become 
the absolute dominate feature in the spectrum. Furthermore, the intensities of different 
spectral features follow the order of v{C-C) (132 mV) < t(CH 2 ) (172 mV) < 5(CH 2 ) (382 
mV) , completely different from the results of the Trl and the Tr2 configurations. The Au-S 
bonding structure of the Chi configuration differs from that of the Trl and the Tr2, resulting 
a noticeable difference in their spectral profiles for the mode z^(Au-S) around 40 mV. The 
changes in molecular conformations seems to be the major cause for the large difference 
in the spectral intensity distributions of two devices. It is found that the molecule in the 
Chi configuration is twisted around the otherwise linearly oriented molecular back-bone in 
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the Trl configuration. It is interesting to note that the spectrum of the Chi configuration 
resembles the experimental IETS of an alkanemonothiol molecule, HS(CH2)sH (CI 1)0, 
quite well, as clearly demonstrated in Figure 0] The molecule-metal bonding structure of 
the Cll is very different from that of octanedithiol. Such a difference should be reflected by 
the spectral profiles of the Au-S modes at the low energy region. The large difference in the 
experimental spectral intensity distribution related to the molecular vibration modes 
implies that the molecular conformations in two experimental setups is very different. The 
molecular back-bone of the octanedithiolate junction in the device of Wang et al.|2f should 
be linear, while it is slightly twisted for the Cll in the device of Kushmerick et al|l|. 

In conclusion, we have proposed a new computational scheme that is capable of describing 
the IETS of molecular junctions with unprecedented accuracy. Our first-principles calcula- 
tions provide reliable assignments for the experimental spectra and reveal important details 
that are not accessible in the experiment. We have also demonstrated that the IETS is a 
powerful characterization tool for molecular devices. 
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TABLE I: Assignments of the vibrational modes observed in the IETS of the octanedithiolate 
junction. The calculated peak positions (in meV), the full width at half maxmimum (FWHM) of 
a single spectral line (in meV), the beginning and the ending of a spectral band (in cm -1 ), and 
the total width of the band A (in meV) are given. The experimental results of Wang et al. |2j are 
also listed for comparison. 
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FIG. 1: Structures of the gold-octanedithiol-gold junctions with triangle (A) and chain (B) local 
contacts. 
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FIG. 2: Inelastic electron tunneling spectrum of the octanedithiol junction from (A) experiment [jj, 

(B) calculation for the Trl configuration and (C) calculation for the Tr2 configuration. The I-V 

curves are given in the inserts of (A) and (B). The working temperature is 4.2 K. 
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FIG. 3: Temperature dependent IETS of the octanedithiolate junction from (A) experiment 
and (B) calculation for the Trl configuration. The intensity is in arbitrary unit. 
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FIG. 4: Inelastic electron tunneling spectroscopy from (A) experiment for Cll monothiol jlj and 
(B) calculation for the Chi configuration. The working temperature is 4.2 K. 
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